Detection of the Power Spectrum of Cosmic Microwave Background Lensing by the 

Atacama Cosmology Telescope 
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We report the first detection of the gravitational lensing of the cosmic microwave background 
through a measurement of the four-point correlation function in the temperature maps made by 
the Atacama Cosmology Telescope. We verify our detection by calculating the levels of potential 
contaminants and performing a number of null tests. The resulting convergence power spectrum at 
2-degree angular scales measures the amplitude of matter density fluctuations on comoving length 
scales of around 100 Mpc at redshifts around 0.5 to 3. The measured amplitude of the signal 
agrees with Lambda Cold Dark Matter cosmology predictions. Since the amplitude of the conver- 
gence power spectrum scales as the square of the amplitude of the density fluctuations, the 4-sigma 
detection of the lensing signal measures the amplitude of density fluctuations to 12%. 



Introduction. — The large-scale distribution of matter 
deflects the paths of microwave background photons by 
roughly 3' [1], a scale larger than the < 1.4' angular reso- 
lution of the Atacama Cosmology Telescope (ACT) . This 
gravitational lensing imprints a distinctive non-Gaussian 
signature on the temperature pattern of the microwave 
sky [2]. Since the cosmic microwave background (CMB) 
temperature fluctuations are very nearly Gaussian [3] 
with a power spectrum now well characterized by WMAP 
[4] and ground-based experiments [5, 6], measurements of 
the distinctive four-point correlation function due to lens- 
ing yield a direct determination of the integrated mass 
fluctuations along the line of sight [2] . 

Previous analyses have detected the lensing signature 
on the microwave sky through cross-correlations of large- 
scale structure tracers with WMAP data [7, 8], or seen 



the signature of lensing in the temperature power spec- 
trum at < 3 (T [5, 6]. Here, we report the first measure- 
ment of the lensing signature using only the CMB tem- 
perature four-point function and constrain the amplitude 
of the projected matter power spectrum. 

Data. — ACT is a six-meter telescope operating in the 
Atacama Desert of Chile at an altitude of 5200 meters. 
The telescope has three 1024-element arrays of super- 
conducting transition-edge sensing bolometers, one each 
operating at 148 GHz, 218 GHz, and 277 GHz. Pre- 
vious ACT team publications describe the instrument, 
observations, and data reduction and initial scientific re- 
sults [6, 9-11]. The analysis presented here is done on a 
324-square-degree stripe of average noise level ~ 23 /iK- 
arcmin, made from three seasons of 148 GHz observations 
of the celestial equator. The region is cut into six equally 
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sized (3x18 degree) patches on which we perform Icnsing 
reconstruction separately, and then combine the results 
with inverse variance weighting. 

The ACT temperature maps (made as in [6]) are fur- 
ther processed to minimize the effects of atmospheric 
noise and point sources. Temperature modes below 
i = 500 as well as a 'stripe' of width i = 180 along the 
Fourier axis corresponding to map declination are filtered 
out to reduce the effects of non-white atmospheric noise 
and scan-synchronous noise respectively [6]. Resolved 
point sources with a signal-to-noise (S/N) greater than 5 
are identified in a match-filtered map [11]. An ACT beam 
template scaled to the peak brightness of each of these 
sources is subtracted from the raw data. Using an algo- 
rithm inspired by the CLEAN algorithm [12], we repeat 
this filtering, source identification, and subtraction un- 
til there are no S/N> 5 identifications. Because the 148 
GHz data also contains temperature decrements from the 
thermal Sunyaev-Zel'dovich (SZ) effect in galaxy clusters, 
the entire subtraction algorithm is also run on the neg- 
ative of the map. The effect of unresolved point sources 
is minimized by filtering out all data above I = 2300. 

Methods. — Gravitational lensing remaps the CMB 
temperature fluctuations on the sky: T{h) = 
T(n -f Q:(n)), where a(n) is the deflection field and un- 
lensed quantities arc denoted by a tilde. In this paper, 
we compute the power spectrum of the convergence fleld, 
K = • a, using an optimal quadratic estimator [13]: 
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where L,L' are coordinates in Fourier space (using 
the flat-sky approximation), g defines filters that can be 
tuned to optimize signal-to-noisc, is a normalization, 
and the second term is the Gaussian part of the four-point 
function. We will refer to the second term as the "Gaus- 
sian bias" , as it is a Gaussian term one subtracts from 
the full four-point function to obtain the non-Gaussian 
lensing signal. We normalize the estimator applying the 
standard formula in [13] using the mean cross-power spec- 
trum estimated from season-splits of the data. 

While the optimal quadratic estimator has the advan- 
tage of maximizing the signal-to-noise, an experimental 
measurement of its amplitude involves subtracting two 
large numbers (the full four-point function and the bias). 
Depending on the quality of data and the relevant length 
scales, this Gaussian four-point bias term can be up to an 
order of magnitude larger than the lensing convergence 
spectrum. As the size of the Gaussian bias term depends 
sensitively on the CMB temperature power spectrum, 
foregrounds and noise, calculating it to sufficient accu- 



racy using the standard simulation or theory approach is 
very difficult, and can lead to large discrepancies. Smidt 
et al. [14] use this standard approach for an analysis 
of the WMAP data, and report a detection significance 
larger than expected from Fisher information theory. An 
alternative approach that does not require this subtrac- 
tion is presented in [15]. 

In this paper, we use the data themselves to obtain a 
first approximation to the Gaussian bias part of the four- 
point function, then compute a small correction using 
Monte Carlo simulations. We first generate multiple ran- 
domized versions of the original data map. The Fourier 
modes of these randomized maps have the same ampli- 
tude as the original map, but with their phases random- 
ized. This destroys any non-Gaussian lensing correlation 
between modes, yet approximately preserves the Gaus- 
sian part of the four point function we wish to model. By 
then averaging the Gaussian biases calculated for many 
realizations of randomized maps, we obtain a good esti- 
mate of the second term in Eq. (1). The small correction 
we subtract from our estimator (a "null bias" at high ^ 
due to spatially varying noise and window functions) is 
easily calculated from Monte- Carlo simulations. A simi- 
lar approach has been suggested by [16]. 

Simulations. — We test our lensing reconstruction 
pipeline by analyzing a large number of simulated lensed 
and unlensed maps. The simulated maps are obtained 
by generating Gaussian random fields with the best 
fit WMAP-I-ACT temperature power spectrum [6, 10], 
which includes foreground models, on maps with the 
ACT pixelization. We then generate lensed maps from 
these unlensed maps by oversampling the unlensed map 
to five times finer resolution, and displacing the pixels 
according to Gaussian random deficction fields realized 
from an input theory. Finally, we convolve the maps 
with the ACT beam, and add simulated noise with the 
same statistical properties as the ACT data, seeded by 
season-split difference maps [6]. 

We apply our lensing estimator to 480 simulations of 
the equatorial ACT temperature map. For each simu- 
lated map we estimate the full four-point function and 
subtract the Gaussian and null bias terms obtained from 
15 realizations of the random phase maps. With 15 real- 
izations, the error on the bias contributes ^ 15% to the 
total error bars. We thus obtain a mean reconstructed 
lensing power spectrum, Eq. (1), as well as the standard 
error on each reconstructed point of the power spectrum. 
The red points in Fig. 1 show the estimated mean con- 
vergence power spectrum from the lensed simulations; it 
can be seen that the input (theory) convergence power 
spectrum is reconstructed accurately by our pipeline. 

Results. — Fig. 2 shows the lensing convergence power 
spectrum estimated from the ACT equatorial data, using 
the value of the Gaussian term as well as the null bias 
obtained from the Monte Carlo simulations previously 
described. The error bars are obtained from the scatter 
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TABLE I. Reconstructed Cf*^ values. 



FIG. 1. Mean convergence power spectrum (red points) from 
480 simulated lensed maps with noise similar to our data. The 
solid line is the input lensing power spectrum, taken from 
the best-fit WMAP+ACT cosmological model. Error bars 
correspond to the scatter of power spectrum values obtained 
from individual maps. 
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FIG. 2. Convergence power spectrum (red points) measured 
from ACT equatorial sky patches. The solid line is the power 
spectrum from the best-fit WMAP-f ACT cosmological model 
with amplitude Al = 1, which is consistent with the measured 
points. The error bars are from the Monte Carlo simulation 
results displayed in Fig. 1. The best-fit lensing power spec- 
trum amplitude to our data is Al ~ 1.16 ± 0.29 



of simulations shown in Fig. 1. 

Here, we introduce the parameter as a lensing con- 
vergence power spectrum amplitude, defined such that 
Al = 1 corresponds to the best-fit WMAP-hACT ACDM 
model (with erg ~ 0.813). The reconstructed points are 
consistent with the theoretical expectation for the con- 
vergence power spectrum. From our results we obtain 
a value of A^ = 1.16 ± 0.29, a 4-cr detection. If we 
restrict our analysis to the first three points, we find 
Al = 0.96 ± 0.31. Fitting our five points to the theory, 
we calculate x^/dof = 6.4/4. Since the lensing kernel 
has a broad peak at z ~ 2 and a conformal distance of 
~ 5000 Mpc, our 4-tT detection is a direct measurement 
of the amplitude of matter fluctuations at a comoving 



FIG. 3. Convergence power spectrum for simulated thermal 
and kinematic SZ maps and point source maps [17] which 
are a good fit to the ACT data. Note that we only show 
the non-Gaussian contribution, as the Gaussian part which 
is of similar negligible size is automatically included in the 
subtracted bias generated by phase randomization. The solid 
line is the convergence power spectrum due to lensing in the 
best-fit WMAP-f ACT cosmological model. 



wavenumber k ~ 0.02Mpc~ around this redshift. 

We estimate potential contamination by point sources 
and SZ clusters by running our reconstruction pipeline on 
simulated patches which contain only IR point sources or 
only thermal or kinetic SZ signal [17], while keeping the 
filters and the normalization the same as for the data 
run. Fig. 3 shows that the estimated spurious conver- 
gence power is at least two orders of magnitude below 
the predicted signal, due partially to our use of only 
temperature modes with £ < 2300. We have also ver- 
ified that reconstruction on simulated maps containing 
all foregrounds (unresolved point sources and SZ) and 
lensed CMB was unbiased. We found no evidence of ar- 
tifacts in the reconstructed convergence power maps. 

Null Tests. — We compute a mean cross-correlation 
power of convergence maps reconstructed from neighbor- 
ing patches of the data map, which is expected to be 
zero as these patches should be uncorrclated. We find 
a x^/dof = 5.8/4 for a fit to zero signal (Fig. 4, upper 
panel). For the second null test we construct a noise 
map for each sky patch by taking the difference of maps 
made from the first half and second half of the season's 
data, and run our lensing estimator. Fig. 4, lower panel, 
shows the mean reconstructed convergence power spec- 
trum for these noise-only maps. Fitting to null we calcu- 
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FIG. 4. Upper panel: Mean cross-correlation power spectrum 
of convergence fields reconstructed from different sky patches. 
The result is consistent with null, as expected. Lower panel: 
Mean convergence power spectrum of noise maps constructed 
from the difference of half-season patches, which is consistent 
with a null signal. The error bars in either case are determined 
from Monte Carlo simulations, and those in the lower panel 
are much smaller as they do not contain cosmic variance. 

late = 5.7 for 4 degrees of freedom. The null test is 
consistent with zero, showing that the contamination of 
our lensing reconstruction by noise is minimal. We also 
tested our phase randomization scheme by randomizing 
the phases on a map, using it to reconstruct a conver- 
gence map, and cross correlating it with a reconstruction 
from the same map but with a different phase random- 
ization; our results were consistent with null as expected. 

Implications and Conclusions. — We have reported a 
first detection of the convergence power spectrum of the 
cosmic microwave background due to gravitational lens- 
ing. The inferred amplitude of the lensing signal is con- 
sistent with theoretical expectations of the basic cosmo- 
logical model. A detection is also anticipated from the 
South Pole Telescope team. Data from the Planck satel- 
lite [18], and CMB polarization measurements with ACT- 
Pol, SPTPol, PolarBear and other next generation exper- 
iments [19] will yield even more accurate measurements of 
CMB lensing. Such measurements are also an important 
goal for a future polarization satellite mission [20]. This 
work is the first step of an exciting research program. 
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